Long-range order in quasi-one-dimensional conductors 
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Abstract 

We study formation of the charge-density wave long-range order in a system of repulsive ID elec- 
trons coupled to 3D phonons. We show that the CDW can be stabilized by interaction with phonons 
in quasi-lD crystals and semiconducting nanowires. In the case of metallic atomic chains, interac- 
tion with phonons of 3D substrate is not enough, and violation of the translational invariance by 
commensurable perturbation or disorder is needed. Possibility of stabilization of superconductivity 
in ID electrons with attraction by means of tunnel coupling to a 3D metal is considered. 

PACS numbers: 71.10.Pm, 71.10.Hf, 71.27.+a, 71.45.Lr, 73.22.-f, 73.63.Nm 



1 



The basic electronic properties of three-dimensional (3D) solids are usually well described 
within Landau's Fermi-liquid picture, including phase transitions into symmetry-broken 
states like the superconducting (SC) or charge- density wave (CDW) state. Even in cases 
where the bare interaction is not weak, energy and momentum conservation reduce the rel- 
evant part of the phase space making a perturbative approach possible. This is not the 
case in ID systems where interaction is always strong. Single-electron quasi-particles do not 
exist, the only low energy excitations are charge and spin collective modes. Such a state is 
called the Luttinger liquid (LL) (for a review see Ref. [1]). A peculiarity of ID systems is 
that long-range order (LRO) cannot exist [2]. In practice however ID systems are embedded 
in a 3D world which results in the coupling of ID electrons to substrates, gates, etc.. Such 
systems are quasi-lD rather than purely ID systems. Then LRO and phase transitions into 
symmetry-breaking states are not excluded. 

It was shown recently [3] that interaction of ID electrons with 3D metallic gate can lead to 
quantum phase transitions and to a power-law decay of the order parameter at zero temper- 
ature. We will show below that coupling to 3D environment can result in phase transitions 
and true LRO even at finite temperatures. Different realizations of the ID electrons embed- 
ded in 3D environment will be considered. The first one is a semiconductor-based quantum 
wire in which dimensionality of the conduction electrons is reduced by dimensional quan- 
tization. Such wires typically have a diameter of the order of 10 nm, hence, they contain 
thousands of primitive cells in the cross-section perpendicular to the wire [4]. The phonon 
system in such wires can be considered to be effectively 3D. Then electron-phonon coupling 
can result in the CDW LRO, and we will determine the conditions under which it will oc- 
cur. Such CDW state is different from the states characterized by quasi-LRO close to a ID 
Wigner crystal [5, 6] expected in ID quantum wires with long-range Coulomb interaction. 

Another interesting class of quasi- ID conductors are highly anisotropic 3D crystals with 
chain-like structure. These materials exhibit commonly metallic behavior and transitions 
to symmetry broken states described in terms of Fermi-liquid ideas. For instance, inorganic 
quasi-lD metals like blue bronzes, NbSe 3 etc. undergo Peierls' transitions to a CDW state. 
LL in such metals is not formed because of the instability towards Fermi-liquid behavior 
around the Fermi energy in the presence of small inter-chain hopping integral t± <C Ep [7]. 
In the same time inter-chain coupling makes them three-dimensional and provides possibility 
of LRO. Whereas the main properties of quasi-lD conductors at higher temperatures are 
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well understood, at low temperatures they demonstrate many intriguing properties which 
are not yet explained convincingly, e.g., the anomalous behavior of the dielectric function 
interpreted as a new glassy phase [8]. Furthermore a transition to a low temperature state 
characterized by LL like conductivity was detected in focused-ion beam processed crystals 
[9]. In order to account for such behavior the possibility of stabilization of the LL state 
by defects in quasi-lD metals was considered [10], but electron-phonon coupling and the 
formation of CDWs was not taken into account. Though formation of LL in quasi- ID metals 
is problematic because of the instability towards Fermi-liquid behavior, LL-like behavior is 
expected to be seen if t± is small compared to the other energy scales of the system related, 
e.g., to disorder or to energy gap in the spectrum. The formation of a CDW energy gap can 
indeed result from the coupling of the quasi-lD crystals electrons to 3D phonons, as will be 
studied below. We also consider a strictly ID electronic system interacting with 3D phonons 
of a substrate (e. g., chain of metallic atoms or a single- wall nanotube). We find that in 
this case the violation of the translational invariance is needed to stabilize the CDW LRO. 
Such a violation can be induced by commensurability effects or disorder in the substrate. 

So far we considered repulsive interaction between electrons. In the opposite case of 
attractive interaction we find that in quasi-lD systems like ID semiconducting nanowires or 
quasi- ID crystals the SC LRO can be stabilized by coupling to a normal 3D metal, while 
in case of strictly ID electronic system this can be performed by a tunnel contact with a 
superconductor. 

For brevity we consider first repulsive spinless electrons exhibiting a ID spectrum de- 
scribed by the Tomonaga-Luttinger model. Further below we will include the modifications 
due to electron spin and the long-range nature of the Coulomb interaction. The electrons 
are assumed to be coupled to 3D phonons. The system is then described by an action 
S = S e \ + + •S'gj.ph consisting of the electronic, the free-phonon part and the electron- 
phonon interaction, respectively. S e \ can be expressed in bosonized form in terms of the 
displacement field & p (t, x, n), where x and r denote space and imaginary time, and n the 
chain index, respectively [11]: 



Here v p = vf/K p is the velocity of plasmons, vf the Fermi velocity, and K p the LL parameter 
measuring the strength of the interaction: K p < 1 for repulsive and K p > 1 for attractive 
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interaction. Sum over n represents summation over conducting chains of quasi-lD crystal 
or over elementary cells for semiconductor quantum wire, say, in Wannier representation. 
Free phonons are described by the elastic displacement field (p with a harmonic action 



Sph — ~ 



S f d 3 q 



dr 



(M 2 + ^(qV 5 



(2) 
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Here S is the cross section area of the primitive cell in the perpendicular plane, oj 2 h (c\) = 
c \q\ + C ±Q±- We need to keep only slowly varying part of the 2k F component of the dis- 
placement field leading to the CDW transition, therefore, in most of the cases the weak 
dependence of oo p h(q) on q\\ can be neglected. Finally, the electron-phonon part reads 



"'el—ph 



= — / dxdr 



(3) 



where 7 is the electron-phonon coupling constant. 

A similar action was used by Voit and Schulz [12] in the study of electron-phonon interac- 
tion in ID conductors. They found modifications of the electron spectrum and the phonon 
dynamics, but no true LRO was present in the purely ID case considered by them. To solve 
the problem for the case of 3D phonons we follow Voit and Schulz and integrate out phonon 
degrees of freedom. This leaves the effective action S e jf = Sgj + with 



>int 



- X! J dxdrdriD^r — Ti,n — rii) 

cos[2$(r, x, n) - 2<E>(tl, x, ni)\ 



(4) 



where D v is the free phonon propagator, its Fourier transform reads D^cu, q_i_) = {oj 2 + 

^(q±)) _1 - 

Next we consider the diagrammatic series for correlation function D p {r, x, n) = 
($ p (r, x, n)$ p (0, 0, 0)). It can be found in a standard way in terms of free correlation func- 
tion £)(°) and mass operator S. 



We calculate S in the lowest approximation in electron-phonon coupling, and find 
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p (0,0)-D p (ri,ni)> 



It is not simple to solve equation (6) in general case, but we can solve it in an interesting 
limiting case when temperature and phonon frequencies are small compared to electronic 
energies. Note that if we neglect fluctuations of the displacement field (p then the action 
Eq. (1,3) reduces to that of the sine-Gordon model. For K p < 1 the latter describes the 
state with the excitations spectrum characterized by a finite dynamic mass M [1]. So we 
expect that the mass operator must yield a mass M in the denominator of the correlation 
function D p (5) at frequencies u ~ M ^> u 2 k F = C||(2fcf). On the other hand, according 
to the form of Eq. (3) the mass operator must vanish in the limit u, q± — > 0. So we seek 
T*(u,q±) in the form 

S(U; ' q±) ~ ~%v^ (uo 2 + c\q\ + CM 2 ) ' (7) 
and assume that the parameters to be found, ( and M, satisfy conditions ( oc (cj 2feF /M) 2 <C 
1. Combining Eqs. (7) and (5) we obtain the correlation function 

itv f (uj 2 + c\q\ + CM 2 ^ 
(u 2 + u 2 )(u 2 + u; 2 ) 

The poles of the denominator describe two eigenmodes instead of a single mode only in 
Ref. [12], the number of modes corresponding to number of variables, <E> P and <p, 

(■ yr 2 „2 2 

2 ~ A/f 2 4- « 2 7, 2 ^ 2 ~ r 2 « 2 4- S ll! p (Q\ 

V ~M +q\\V p , u 2 ~ c ± g ± + M 2 + g 2^2 • W 



~ I, a , , ,2\/, ,2 i ,.2\ • W 



The mode a>i is related mainly to electronic fluctuations, while the soft mode u 2 mainly 
controls phase fluctuations of the new Peierls-deformed lattice, i.e. of the CDW. Note that 
the infrared divergence in (<3> 2 ) = / ^^^D p (uj, q) is removed now and the expectation 
values of the 2/c^-component of the electronic density is finite 

o 1 ( M \ Kp (m) 

So using Eq. (8) in (6) we find that the mass operator has indeed the form (6) provided 
K p is not too close to 1 (i.e. for strong enough electron-electron repulsion) and 

7> 2feF «M, In— » VCln 11 ™ ■ (11) 



with values of parameters 



/ /o\V \ V(i-^p) , ,2 



where A = — «C 1 is dimensionless electro n-phonon coupling constant, and A ps 1/a is 
a large momentum cut-off. The last relation in Eq. (11) is easily satisfied if the anisotropy of 
the phonon spectrum is not too large, but cannot be satisfied in case of ID phonon spectrum. 

When this result is applied to quasi-lD crystals there is another mechanism of three- 
dimensionality, namely, inter-chain tunnelling that could lead to instability of the LL be- 
havior. However, in the presence of finite mass M induced by coupling to phonons contribu- 
tion of the inter-chain hopping can be neglected provided the inter-chain coupling is small 
enough. 

The same results can be obtained by neglecting fluctuations of the amplitude of the 2k F - 
displacement field, ip = r]e lx , and using the self-consistent harmonic approximation (SCHA) 
[1]. In this way one can calculate (x 2 ) 



(X 



2\ 



2^K P J ff$^K ( ^ ) at T< u ph , 
ttK p V( f S ( fU Tu!2kF atT» Wph . 

^ (^Y c ± q ± uj ph P a 



where the elliptic integral K « l n (^hE_\ f or q± _^ q Thus CDW phase fluctuations, 

(x 2 ), are small under conditions (11). If we neglect small value (x 2 ) then our problem is 
reduced to the exactly solved sine-Gordon and massive Tirring models [1, 12]. The value 
of M (12) agrees with these models deep in the massive region (when interaction is strong 
enough). For non- interacting electrons, K — > 1, the exact solutions yield the result of the 
classical Peierls theory M ~ Af^e^x. 

A finite value of the order parameter in Eq. (10) follows from dependence of the mass 
operator (7) on q±. In a strictly ID electronic system interacting with 3D phonons of a 
substrate the dependence on q± is absent and S = at u>, q± — > 0. This is because the phase 
fields $ or x themselves do not couple to physical fields in a translationally invariant system. 
Indeed, in the limit u, q — > they describe translations of the electrons and the lattice, 
respectively, as a whole. The physical fields are related to time and spatial derivatives of the 
phase fields. Hence phase fluctuations cannot diffuse into 3D space effectively because their 
interaction with 3D phonons disappears in the infrared limit. Therefore the contribution 
of the mode 002 to (x 2 ) remains divergent and coupling to 3D phonons does not lead to 
LRO similar to results of Refs. [12]. To remove the infrared divergence in (x 2 ) one needs to 
break the translational invariance, e.g., by a commensurate potential in the substrate or by 
disorder. 
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Let us first consider the violation of translational invariance by a potential with period 
7t//cf along the x-direction. The additional contribution to the action, 

5S = TJ2[ dxdT[0* 2kF <p(n) + 2kF <p*(n)], (14) 

n J 

is presented in a form resembling the electron-phonon interaction (3). Here 2 k F describes 
density perturbations, 2kF = \0 2 k F \ exp(i2kpx + ixo), inducing a potential, T0 2kF , that 
acts on 2/ci?-phonons. 

Taking into account that destruction of the LRO is controlled by phase fluctuations one 
can neglect fluctuations of the amplitude of the phonon field [12]. Then using the SCHA 
we integrate out the electronic degrees of freedom and the phonon operators at all chains 
except the metallic one, ip(n = 0) = r\e %x . This results in an effective action for the CDW 
phase 

SS eff = T J dxdr\0 2kF \r] cos(x - Xo)- (15) 

With Eq. (15) the results for single metallic chain become similar to (12-13) if we substi- 
tute c\q\rf by T\0 2 k F \i]. Then instead of Eq. (13) we obtain 

In the limit of small commensurate potential, T <C 7 or 2 k F \fS <C 1 the logarithm can be 
large. However, even in this case (x 2 ) can be small due to the small factor y^. 

The divergence in (x 2 ) can be removed also by disorder. In this case one has to substitute 
2 k F in Eq. (14) by the local impurity potential in the substrate V(r) = J2i V5(r — r,j). This 
leads to a term similar to Eq. (15) but with J2i S(x — Xi) cos(x — 2kpXi) replacing cos(x — Xo)- 
Though this will remove the divergence in (x 2 ) due to quantum fluctuations, true LRO is 
not expected since the latter is destroyed by the disorder. Nevertheless, for weak disorder 
the phase coherence length will be quite large as it is in CDW conductors. 

For system of spinful electrons the results are qualitatively similar to the spinless case. 
The spin is taken into account by adding to Eq.(l) a spin term analogous to the charge 
one [1]. The electron-phonon term (3) must be also modified [12], namely, the integrand is 
multiplied by cos \/2& a where $ CT is the spin phase field, and factors 2 in the exponents are 
substituted by y/2. After going through the same steps as in the spinless case we find mass 
operators for correlation functions D p and D a . For the charge sector S p has a form similar 
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to Eq. (7), while S CT is not dispersive, so the spin excitations are described by dispersion 
relation similar to uj\ (9) but with different velocity. The masses for charge and spin sectors 
are equal in the main approximation 

M = Av F [\ Kf» K«°/(aA) 2 } 1/(2 - K p- K °\ (16) 

Since we consider repulsive electrons we can ignore the backscattering in the interaction. 
While in case of inter-electronic attraction it leads to the spin-gapped Luther-Emery phase 
(c/. [1]), in our problem it results only in small corrections. 

Up to now we considered the short-range interaction and treated K p as a constant. In the 
case of Coulomb interaction it depends on wave vector. The exact form of K p (q) depends on 
crystal structure, but in the long-wavelength limit one can describe interaction by the Fourier 
transform of the Coulomb potential. Then K 2 (q) = (1 + n 2 /q 2 ) where 1/k is Thomas-Fermi 
screening length. This results in qualitative modifications of the spectrum of the eigenmodes. 
Inserting v p = vf/K p into (9) we find that the mode lo\ approaches at q±_ = the plasma 
frequency, and the phason mode u 2 is not a soft mode any more. Similar hardening occurs 
in the LL density mode without the coupling to phonons [11] and in the CDW phase mode 
in classical Peierls theory [13]. 

Similar to the CDW case, but with action written in terms of dual field 6 P [1], one can 
study stabilization of SC LRO in a ID electronic system with attraction, (K p > 1). Since 
the SC phase is conjugate to the particle number operator and the fixed number of particles 
makes fluctuations of the phase infinite one can try to let the particle number to fluctuate 
in order to diminish fluctuations of the SC phase by coupling to a normal or SC 3D metal 
via a tunnel junction providing a way to transfer fluctuations of the phase to 3D space. We 
describe the coupling to the 3D metal by a tunnel Hamiltonian. The 3D metal is presented by 
the BCS Hamiltonian in the Bogolyubov self-consistent approximation. Then we integrate 
out the fermions of the 3D metal and obtain an effective action describing coupling to the 
SC order parameter, A, of the 3D metal: 

Stl f) = E / dxdrt(e 2 ^A + e^ e A*). (17) 

We describe A by a Ginzburg-Landau free energy for a 3D superconductor with critical 
temperature T c . So we consider coupling to a superconductor for T < T c , and to a normal 
metal for T > T c . We find that the SC LRO can be stabilized by coupling to a normal 3D 
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metal in quasi-lD structures where dependence of P on transverse coordinates is possible. 
In strictly ID systems coupling to a normal metal does not stabilize the LRO, and only 
quantum phase transition as discussed in Ref. [3] are possible. This happens because, like 
in the CDW case, physical fields are related to derivatives of the SC phase, and coupling 
to 3D space vanishes in the infrared limit. So, again, the symmetry breaking is needed. 
The symmetry is violated if the 3D metal is in the SC state, (A) 7^ in Eq. (17), and the 
problem is reduced to the sine-Gordon model. Hence, the LRO becomes possible. 

In conclusion, we found that coupling of electrons in quasi-lD systems to 3D objects 
should lead to symmetry-breaking transitions described in terms of the LL ideas. The 
transition to the CDW state was studied in case when phonon frequency is much smaller 
than typical electronic energy and the repulsion between electrons is strong enough. Both 
conditions are satisfied in quasi-lD crystals where typical values of the Peierls gap are above 
500K, and the interaction parameter K p is quite small, K p ~ 0.1 4- 0.3 [10]. Our study of 
the SC transition explains possibility of the LRO formation in quasi-lD systems but it has 
more model character because we consider electronic attraction as granted and do not study 
mechanism of attraction. 
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Research. 
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